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Abstract 

By using a simple analytical model based on counting random multiple impulses inflicted 
on photons by a network of cosmic strings we show how to construct the general q-point tem- 
perature correlation function of the Cosmic Microwave Background radiation. Our analysis 
is sensible specially for large angular scales where the Kaiser- Stebbins effect is dominant. 
Then we concentrate our study on the four-point function and in particular on its zero-lag 
limit, namely, the excess kurtosis parameter, for which we obtain a predicted value of ~ 10~ 2 . 
In addition, we estimate the cosmic variance for the kurtosis due to a Gaussian fluctuation 
field, showing its dependence on the primordial spectral index of density fluctuations n and 
finding agreement with previous published results for the particular case of a flat Harrison- 
Zel'dovich spectrum. Our value for the kurtosis compares well with previous analyses but 
falls below the threshold imposed by the cosmic variance when commonly accepted param- 
eters from string simulations are considered. In particular the non-Gaussian signal is found 
to be inversely proportional to the scaling number of defects, as could be expected by the 
central limit theorem. 
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1 Introduction 



A central concept for particle physics theories attempting to unify the fundamental interac- 
tions is the concept of symmetry breaking. This symmetry breaking plays a crucial role in 
the Weinberg-Salam standard electroweak model ( Masiero 1984) whose extraordinary suc- 
cess in explaining electroweak scale physics reaches a precision rarely found before in other 
areas of science (Koratzinos 1994). In the context of the standard Big Bang cosmological 
theory the spontaneous breaking of fundamental symmetries is realized as a phase transition 
in the early universe. Such phase transitions have several exciting cosmological consequences 
thus providing an important link between particle physics and cosmology. 

A particularly interesting cosmological issue is the origin of structure in the universe. 
This structure is believed to have emerged by the gravitational growth of primordial matter 
fluctuations which are superposed on the smooth background required by the cosmological 
principle, the main assumption of the Big Bang theory. 

The above mentioned link of cosmology to particle physics theories has led to the gener- 
ation of two classes of theories which attempt to provide physically motivated solutions to 
the problem of the origin of structure in the universe. According to the one class of theories, 
based on inflation, primordial fluctuations arose from zero point quantum fluctuations of a 
scalar field during an epoch of superluminal expansion of the scale factor of the universe. 
These fluctuations may be shown to obey Gaussian statistics to a very high degree and to 
have an approximately scale invariant power spectrum. 

According to the second class of theories, those based on topological defects, primordial 
fluctuations were produced by a superposition of seeds made of localized energy density 
trapped during a symmetry breaking phase transition in the early universe. Topological 
defects with linear geometry are known as cosmic strings and may be shown to be consis- 
tent with standard cosmology unlike their pointlike (monopoles) and planar (domain walls) 
counterparts which require dilution by inflation to avoid overdosing the universe. Cosmic 
strings are predicted to form during a phase transition in the early universe by many but 
not all Grand Unified Theories (GUTs). 

The main elegant feature of the cosmic string model that has caused significant attention 
during the past decade is that the only free parameter of the model (the effective mass 
per unit length of the wiggly string /j, ) is fixed to approximately the same value from two 
completely independent directions. From the microphysical point of view the constraint 
Gfi ~ (jtigut l m p) 2 — 10~ 6 is imposed in order that strings form during the physically 
realizable GUT phase transition. From the macrophysical point of view the same constraint 
arises by demanding that the string model be consistent with measurements of Cosmic 
Microwave Background (CMB) anisotropies and that fluctuations are strong enough for 
structures to form by the present time. This mass-scale ratio for Gfi is actually a very 
attractive feature and even models of inflation have been proposed (Freese et al. 1990, Dvali 
et al. 1994) where observational predictions are related to similar mass scale relations. 

Cosmic strings can account for the formation of large scale filaments and sheets (Vachas- 
pati 1986; Stebbins et al. 1987; Perivolaropoulos, Brandenberger & Stebbins 1990; Vachas- 
pati & Vilenkin 1991; Vollick 1992; Hara & Miyoshi 1993), galaxy formation at epochs 
z ~ 2 — 3 (Brandenberger et al. 1987) and galactic magnetic fields (Vachaspati 1992b). 
They also generate peculiar velocities on large scales (Vachaspati 1992a; Perivolaropoulos & 
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Vachaspati 1993), and are consistent with the amplitude, spectral index (Bouchet, Bennett & 
Stebbins 1988; Bennett, Stebbins & Bouchet 1992; Perivolaropoulos 1993a; Hindmarsh 1993) 
and the statistics (Gott et al. 1990; Perivolaropoulos 1993b; Moessner, Perivolaropoulos & 
Brandenberger 1993; Coulson et al. 1993; Luo 1994) of the CMB anisotropies measured by 
COBE on angular scales of order 6 ~ 10°. 

Strings may also leave their imprint on the CMB mainly in three different ways. The 
best studied mechanism for producing temperature fluctuations on the CMB by cosmic 
strings is the Kaiser-Stebbins effect (Kaiser & Stebbins 1984; Gott 1985). According to this 
effect, moving long strings present between the time of recombination t rec and today produce 
(due to their deficit angle (Vilenkin 1981)) discontinuities in the CMB temperature between 
photons reaching the observer through opposite sides of the string. Another mechanism for 
producing CMB fluctuations by cosmic strings is based on potential fluctuations on the last 
scattering surface (LSS). Long strings and loops present between the time of equal matter and 
radiation t eq and the time of recombination t rec induce density and velocity fluctuations to 
their surrounding matter. These fluctuations grow gravitationally and at t rec they produce 
potential fluctuations on the LSS. Temperature fluctuations arise because photons have 
to climb out of a potential with spatially dependent depth. A third mechanism for the 
production of temperature anisotropies is based on the Doppler effect. Moving long strings 
present on the LSS drag the surrounding plasma and produce velocity fields. Thus, photons 
scattered for last time on these perturbed last scatterers suffer temperature fluctuations due 
to the Doppler effect. 

It was recently shown (Perivolaropoulos 1994) how, by superposing the effects of these 
three mechanisms at all times from t rec to today, the power spectrum of the total temperature 
perturbation may be obtained. It turns out that (assuming standard recombination) both 
Doppler and potential fluctuations at the LSS dominate over post-recombination effects on 
angular scales below 2°. However this is not the case for very large scales (where we will 
be focusing in the present paper) and this justifies our neglecting the former two sources of 
CMB anisotropies. The main effect of these neglected perturbations is an increase of the 
gaussian character of the fluctuations on small angular scales. 

The main assumptions of the model were explained in (Perivolaropoulos 1993a). Here we 
will only review them briefly for completeness. As mentioned above, discontinuities in the 
temperature of the photons arise due to the peculiar nature of the spacetime around a long 
string which even though is locally flat, globally has the geometry of a cone with deficit angle 
8nGfj,. Several are the cosmological effects produced by the mere existence of this deficit 
angle (Shellard 1994), e.g., arcsecond-double images from GUT strings at redshifts z ~ 1, 
flatten structures from string wakes or elongated filamentary structures from slowly moving 
long wiggly strings and, of more relevance in our present work, post-recombination CMB 
anisotropy (White, Scott and Silk 1994) string induced effects (Kaiser & Stebbins 1984). 

The magnitude of the discontinuity is proportional not only to the deficit angle but also 
to the string velocity v s and depends on the relative orientation between the unit vector 
along the string s and the unit photon wave-vector k. It is given by (Stebbins 1988) 

AT 

— = ±AnGfiv s -f s k ■ (v s x s) (1) 
where 7 S is the relativistic Lorentz factor and the sign changes when the string is crossed. 
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Also, long strings within each horizon have random velocities, positions and orientations. 

We discretize the time between t rec and today by a set of iV Hubble time-steps ij such 
that £j + i = tid, i.e., the horizon gets multiplied by 5 in each time-step. For z rec ~ 1400 we 
have N ~ log (5 [(1400) 3/2 ]. 

In the frame of the multiple impulse approximation the effect of the string network on 
a photon beam is just the linear superposition of the individual effects, taking into account 
compensation (Traschen et al. 1986; Veeraraghavan & Stebbins 1990; Magueijo 1992), that 
is, only those strings within a horizon distance from the beam inflict perturbations to the 
photons. 

In the following section we show how to construct the general q-point function of CMB 
anisotropies at large angular scales produced through the Kaiser- Stebbins effect. Explicit 
calculations are performed for the four-point function and its zero-lag limit, the kurtosis. 
Next, we calculate the (cosmic) variance for the kurtosis assuming Gaussian statistics for 
arbitrary value of the spectral index and compare it with the string predicted value (section 
3). Finally, in section 4 we briefly discuss our results. 



2 The Four— Point Temperature Correlation Function 

According to the previous description, the total temperature shift in the 7 direction due to 
the Kaiser-Stebbins effect on the microwave photons between the time of recombination and 
today may be written as 

A rp N M 

—(7) = 47^.7. E E r n {i) (2) 

n=l m=l 

where /3 mn (7) gives us information about the velocity v mn and orientation s mn of the mth 
string at the nth Hubble time-step and may be cast as f3 mn {^) = 7 ■ R mn , with R mn = 
v mn x s mn = (sin 9 mn cos(p mn , sin 9 mn sin mn , cos9 mn ) a unit vector whose direction varies 
randomly according to the also random orientations and velocities in the string network. In 
Eq.(^), M denotes the mean number of strings per horizon scale, obtained from simulations 
( Allen & Shellard 1990, Bennett & Bouchet 1988) to be of order M ~ 10. 

Let us now study the correlations in temperature anisotropies. We are focusing here in 
the stringy-perturbation inflicted on photons after the time of recombination on their way 
to us. Photons coming from different directions of the sky share a common history that may 
be long or short depending on their angular separation. We therefore take 9 tp ~ zitp)^ 1 ^ 2 
to be the angular size of the horizon at Hubble time-step t p (1 < p < N). The kicks 
on two different photon beams separated by an angular scale au = arccos(7i greater 
than 9t p will be uncorrelated for the time-step t p (no common history up until ~ t p ) but 
will eventually become correlated (and begin sharing a common past) afterwards when the 
horizon increases, say, at t p > when 9 t , > cti2, i.e., when a%2 fits in the horizon scale. 

Much in the same way, kicks inflicted on three photon beams will be uncorrelated if at 
time t any one of the three angles between any two directions (say, ai 2 , ct23, 0^13) is greater 
than 9t, the size of the horizon at time t. So, in this case, we will be summing (cf. Eq.@) 
over those Hubble time-steps n greater than p, where p is the time-step when the condition 
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9 t = Max[«i2, «23, CK13] is satisfied. The same argument could be extended to any number 
of photon beams (and therefore for the computation of the q-point function of temperature 
anisotropies)^. 

Let us put the above considerations on more quantitative grounds, by writing the mean 
q-point correlation function. Using Eq.(|2|) we may express it as 

AT AT N M 

<—{!!)...—(%)>=? £ E <PT ini -PT qnq > ( 3 ) 

n\,...,n q =l mi,...,m q =l 

where 71. ..j q are unit vectors denoting directions in the sky and where pj 1 ^ 3 = ^yR m i n J with 
7j = (sin 8j cos (f)j, sin 9j sin <f>j, cosOj). In the previous equation we defined £ = AnGnv s ^ s . 
For simplicity we may always choose a coordinate system such that 6\ = 0, 4>i — (i.e., 71 lies 
on the z axis) and 02 = tt/2. The seemingly complicated sum (|3J) is in practice fairly simple 
to calculate because of the large number of terms that vanish due to lack of correlation, after 
the average is taken. The calculation proceeds by first splitting the product (3™ ini ...(3™ qUq 
into all possible sub-products that correspond to correlated kicks (i.e., /3 OTj '" 3 ''s with the same 
pair (rrij,nj)) at each expansion Hubble-step and then evaluating the ensemble average of 
each sub-product by integration over all directions of R. To illustrate this technique we 
evaluate the two, three and four-point functions below. 

The calculation of the two-point function has been performed in (Perivolaropoulos 1993a) 
but we briefly repeat it here for clarity and completeness. Having a correlated pair of 
beams in 71 and 72 directions from a particular time-step p onwards means simply that 
i? mini = R m2n ' 2 <^> n = Hi = n 2 > p, m = mi = ni2] otherwise the i?'s remain uncorrelated 
and there is no contribution to the mean two-point function. Therefore we will have 

AT AT N M 

-^(71)^(72) = e e E (71 • R mn )(i2 ■ R mn ) (4) 

n=p m=l 

where we wrote just the correlated part on an angular scale ai 2 = arccos(7i • 72). The 
uncorrelated part on this scale will vanish when performing the ensemble average (■). Thus 
the mean two-point function may be written as 

AT AT 

< _(7 1 )_(7 3 ) >= e < > N cor {a 12 ) (5) 

where 

N cor (a 12 ) = M[N-3 log 5 (l + ^)] (6) 

tree 

is the number of correlated 'kicks' inflicted on a scale after this scale enters the horizon. 
8 rec is the angular scale of the horizon at the time of recombination. Since we may always 
take 

<3\ = (0, 0, 1) • R /3 2 = (sinai2, 0, COSOJ12) • -R (7) 



2 A general expression (suitable to whatever angular scale and to whatever source of temperature fluc- 
tuations) for the three-point correlation function was given in (Gangui, Lucchin, Matarrese and Mollerach 
1994). Although of much more complexity, similar analysis may be done for the four-point function and a 
completely general expression in terms of the angular trispectrum may be found (Gangui & Perivolaropoulos 
1994). 
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with R = (sin 9 cos 0, sin 9 sin 0, cosQ), < ^(71)^(72) > may be calculated by integrating 
over (9, 0) and dividing by An. The result is 

AT.„ ,A7\ . ~cosai2 AT , N 
< — ( 7l ) — ( 72 ) >= e^^N cor {ai2) (8) 

which may be shown (Perivolaropoulos 1993a) to lead to a slightly tilted scale invariant 
spectrum on large angular scales. 

The three-point function may be obtained in a similar way. However, the superimposed 
kernels of the distribution turn out to be symmetric with respect to positive and negative 
perturbations and therefore no mean value for the three-point correlation function arises (in 
particular also the skewness is zero). On the other hand the four-point function is easily 
found and a mean value for the excess kurtosis parameter (Gangui 1994b) may be predicted, 
as we show below. 

In the case of the four-point function we could find terms where (for a particular time- 
step) the photon beams in directions 71, 72 and 73 are all correlated amongst themselves but 
not the beam 74, which may be taken sufficiently far apart from the other three directions. 
In such a case we have 

,A7\ ,AT, A .AT.. sATV ,A7V AT., ,AT,„ ,A7V ,. 

(—(71)^(72)^(73)^(74)) — ► (—(71)^(72)^(73)) (-jtKIa)) (9) 

and (cf. (Perivolaropoulos 1993a)) this contribution vanishes. 

Another possible configuration we could encounter is the one in which the beams are 
correlated two by two but no correlation exists between the pairs (for one particular time- 
step). This yields three distinct possible outcomes, e.g., ( 4^(7 1 )- A ^(7 2 )) (^(73)^(74)), 
and the other two obvious combinations. 

The last possibility is having all four photon beams fully correlated amongst themselves 
and this yields (^-(71) ^(72)^?- (73)^(74))^ where the subscript c stands for the connected 
part. 

In a way completely analogous to that for the two-point function, the correlated part for 
the combination of four beams gives 

AT AT AT AT N M 

-^(71)^(72)^(73)^(74) = e 4 E T,(^-R mn )(i2-R mn )(i3-R mn )(^-R mn ) (io) 

n=p m=l 

Now we are ready to write the full mean four-point function as 



AT AT AT AT 

' AT „ AT , . N , , AT „ AT n 
(— (7i)^(7 2 ))(^(73)^(74)) + 2terms 



AT AT AT AT 

+( — (71)^(72)^(73)^(74)), 



T T x ' T T 
1 .0 I 2 ,AT_ , AT ^ V AT /A , AT 



3 e 2 iV cor (^) cos(^) 



+ (—(7i)^(72)^(73)^(74)) c (11) 



where for simplicity we wrote this equation using the same scale 9 for all directions on the 
sky (first term after the equality sign); after all, we will be interested in the zero lag limit, 
in which case we will have 9 — > 0. By using Mathematica (Wolfram 1991) it is simple to 
find that the second term includes just the sum of twenty-one combinations of trigonometric 
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functions depending on 9i, 4>i, with i = 1,2,3,4, the spherical angles for the directions % on 
the sky (cf. Eq. ([T0|) ). These terms are the only ones which contribute non-vanishingly after 
the integration over the angles Q mn ^ rnn i s performed (assuming ergodicity). 

When we take the zero lag limit (i.e., aiming for the kurtosis) the above expression gets 
largely simplified. After normalizing by the squared of the variance <r 4 = [^C, 2 N cor (0)] 2 and 
subtracting the disconnected part we get the excess kurtosis parameter 

/C = ^<^T(7i)> - 3 = ^ ~ (1.125 t° llt) x 10- 2 (12) 

where we took M ~ 10 and 5 = 2 implying N ~ 16. Quoted errors take into account possible 
variation of 5 in the range 1.5 < 8 < 3. 

Note that for values of the scaling solution M increasing (large number of seeds) the 
actually small non-Gaussian signal /C gets further depressed, as it could be expected from 
the Central Limit Theorem. 



3 The r.m.s. excess kurtosis of a Gaussian field 

The previous section was devoted to the computation of /C, the excess kurtosis parameter, 
as predicted by a simple analytical model in the framework of the Cosmic String scenario. 
We might ask whether this particular non-Gaussian signal has any chance of actually being 
unveiled by current anisotropy experiments. Needless to say, this could provide a significant 
probe of the structure of the relic radiation and furthermore give us a hint on the possible 
sources of the primordial perturbations that, after non-linear evolution, realize the large scale 
structure of the universe as we currently see it. 

However, as it was realized some time ago (Scaramella & Vittorio 1991; Abbott & Wise 
1984; Srednicki 1993; Gangui 1994a), the mere detection of a non-zero higher order correla- 
tion function (e.g., the four-point function) or its zero lag limit (e.g., the kurtosis) cannot be 
directly interpreted as a signal for intrinsically non-Gaussian perturbations. In order to tell 
whether or not a particular measured value for the kurtosis constitutes a significant evidence 
of a departure from Gaussian statistics, we need to know the amplitude of the non-Gaussian 
pattern produced by a Gaussian perturbation field. Namely, we need to know the r.m.s. 
excess kurtosis of a Gaussian field. 

Let us begin with some basics. Let us denote the kurtosis K = J -^-^f- iff) and assume 
that the underlying statistics is Gaussian, namely, that the multipole coefficients a™ are 
Gaussian distributed. Thus, the ensemble average for the kurtosis is given by the well- 
known formula: (K) = 3cr 4 , where a 2 is the mean two-point function at zero lag, i.e., the 
CMB variance as given by 

a 2 = (C7 2 (0)> = ^E^f^Q 2 ^ 2 (13) 

where the Ci coefficients (here normalized to C^=2 = 1) are also dependent on the value for 
the primordial spectral index of density fluctuations n and are given by the usual expression 
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in terms of Gamma functions (Bond & Efstathiou 1987; Fabbri, Lucchin & Matarrese 1987) 



r(i + f-|)r(f- 
r(€+§-§)r(| + 2) 



(14) 



Q = (Q2)V2 is the rms quadrupole, simply related to the quantity Q rms -ps defined in 
(Smoot et al. 1992; Bennett et al. 1994) by Q = \Z4jrQ rms -ps/T , with mean temperature 
T = 2.726 ± 0.01-K" (Mather et al. 1994). In the previous expression W<? represents the 
window function of the specific experiment. Setting Wo = Wi = automatically accounts 
for both monopole and dipole subtraction. In the particular case of the COBE experimental 
setup we have, for i > 2, W e ~ exp + 1)(3.2 ) 2 ] , where 3.2° is the dispersion of the 

antenna-beam profile, which measures the angular response of the detector (e.g. Wright et 
al. 1992). Sometimes the quadrupole term is also subtracted from the maps and in that case 
we also set W2 = 0. 

However, (K) is just the mean value of the distribution and therefore we cannot know 
its real value but within some error bars. In order to find out how probable it is to get this 
value after a set of experiments (observations) is performed, we need to know the variance 
of the distribution for K. In other words, we ought to know how peaked the distribution is 
around its mean value (K). The width of the distribution is commonly parameterized by 
what is called the cosmic variance of the kurtosis 

o*cv = (K 2 ) - (K) 2 (15) 

It is precisely this quantity what attaches theoretical error bars to the actual value for the 
kurtosis. Therefore, we may heuristically express the effect of a cv on K Gauss as follows: 
K Gauss ~ (K) ± Gov-, at one sigma level (a good approximation in the case of a narrow 
peak). Re-arranging factors we may write this expression in a way convenient for comparing 
it with /C as follows 

v~ K Gauss 1 @CV / t n\ 

ICcv = - 3 ^ ±-^- (16) 

where ICcv is the excess kurtosis parameter (assuming Gaussian statistics) purely due to the 
cosmic variance. Not only is ICcv i n general non-zero, but its magnitude increases with the 
theoretical uncertainty ( acv ) due to the limitation induced by our impossibility of making 
measurements in more than one Universe. 

This gives a fundamental threshold that must be overcome by any measurable kurtosis 
parameter in order for us to be able to distinguish the primordial non-Gaussian signal from 
the theoretical noise in which it is embedded. In other words, unless our predicted value for 
/C exceeds ICcv we w iU n °t be able to tell confidently that any measured value of IC is due 
to primordial non-Gaussianities. 

Now that we know the expression for the cosmic variance of the kurtosis, let us calculate 
it explicitly. We begin by calculating (K 2 ) as follows 

^ = /^/^<fW£\*>> (17) 



8 



By assuming Gaussian statistics for the temperature perturbations ^-(7) we may make use 
of standard combinatoric relations, and get 

AT 4 AT 4 AT 2 . 2/ AT 2 2 

(— (7i)-jT (72)) = 9(— (71)) {-jr (72)) 

AT 2 AT 2 AT AT AT AT 

+ 72(— (7i)><— (7 2 ))(^(7i)^(7 2 )) 2 + 24( — ( 7l ) — ( 72 )) 4 (18) 

As we know, the ensemble averages are rotationally invariant and therefore (^?- 2 (7i)) = 
(C2(0)) = a 2 is independent of the direction 71. Plug this last equation into Eq. (|17"D and we 

get 

(K 2 ) = 9a s + 36a 4 J^d cos a (C 2 (a)) 2 + 12 j dcosa (C 2 (a)) 4 . (19) 

The above integrals may be solved numerically. Then, using this result in the expression for 
a^ v we get JCcv, the value for the excess kurtosis parameter of a Gaussian field. 

It is also instructive to look at Eq . (|19|) in some more detail, so that the actual dependence 
on the spectral index becomes clear. By expanding the mean two-point correlation functions 
within this expression in terms of spherical harmonics and after some long but otherwise 
straightforward algebra we find 



JC, 



cv 



Et(2£ + l)C 2 W 4 
[Ei(2£ + l)C e W 2 } 2 



x 21 1/2 

/TT 4 Y^i r 1/\;2 1 [ST* A^^rn 1 ;m2,m 3 +m in ,rn 3 ,m i ,-m 3 -m i ^ 



[E,(2£+1)C,W ( 
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(20) 



where the coefficients H^£^ = J dtt^Y™ 1 (7)^ {l)Y% 3 \l) , which can be easily ex- 
pressed in terms of Clebsch-Gordan coefficients (Messiah 1976), are only non-zero if the 
indices £ i} rrii (i = 1, 2, 3) fulfill the relations: \£j — £k\ < U < |^j+^fc|> ^1 + ^2 + ^3 = even and 
mi + m 2 + m 3 = 0. In the above equation the n-dependence is hidden inside the multipole 
coefficients C £ (cf. Eq.®). 

Eq.(|20|) shows an analytic expression for computing JCcv which, in turn, represents a 
fundamental threshold for any given non-Gaussian signal. For interesting values of the 
spectral index (say, between 0.8 < 77 < 1.3) we find no important variation in JCcv, being its 
value consistent with the Monte-Carlo simulations performed by Scaramella & Vittorio 1991 
(SV91 hereafter) , see below. These authors concentrated on a Harrison-Zel'dovich spectrum, 
considered also the quadrupole contribution and took a slightly different dispersion width 
for the window function (3.0° in their simulations). 

We coded an IBM RISC 6000 for solving Eq.(^) numerically. We first calculated Kqv for 
a somewhat reduced range for the multipole index £ with quadrupole subtracted, 3 < £ < 5, 
so that we were able to test the result with that obtained by using Mathematica. We found 
perfect agreement and a value fCcv — 1-85 (for n = 1). We may even go further and include 
the quadrupole, i.e. 2 < £ < 5, and in this case we get JCcv — 1.92. This slightly larger 
value obtained after including £ = 2 is not a new feature and simply reflects the intrinsic 
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theoretical uncertainty of the lowest order multipoles; see e.g. (Gangui et al. 1994) for a 
similar situation in the case of the skewness. 

Of course ICcv — 1-85 (or 1.92 with quadrupole) are still a factor 6 above the value ~ 0.3 
found in SV91 (where they included the quadrupole). Our analytical analysis expresses ICcv 
as a ratio of averages (rather than the average of a ratio as in SV91) and therefore exact 
agreement should not be expected. But still the main reason for the discrepancy lies in the 
small range for the multipole index I. We therefore increased the value of £ ma x and checked 
that ICcv monotonically became smaller and smaller. When we were in the range 3 < £ < 10 
we got ICcv — 1-16; instead, for 3 < £ < 15 we got ICcv — 0.98 (for n = 1) -clearly a sensible 
decrease. CPU-time limitations prevent us from carrying out the numerical computations for 
larger values of £ max (usually £ max is chosen of order 30, i.e., beyond those values of £ where 
the exponential suppression of the We makes higher £ contribution to the sums negligible), 
but still we expect ICcv to keep on steadily decreasing (we checked that the decrease in 
ICcv was tinier as £ max got larger). This, together with some previous experience in similar 
computations (Gangui et al. 1994) makes us believe that in the case appropriately large 
values for £ max were used , the SV91 value above mentioned would be attained (taking into 
account the quadrupole subtraction, of course). 

In addition, we checked that a different value of the spectral index n does not change 
the essence of the above considerations. We explored the cosmologically interesting range 
0.8 < n < 1.3 and got values 1.22 > IC CV > 1-08 (taking 3 < £ < 10) and values 1.05 > 
ICcv > 0.88 (taking 3 < £ < 15). We plot ICcv versus the spectral index for £ max = 15 in 
Fig.l. 

Note the small rate of variation of ICcv with n and that, as expected, ICcv takes larger 
values for smaller spectral indexes. This is clearly due to the fact that a small n generates 
more power on large scales (i.e., small £) and precisely these scales are the ones that contribute 
the most to the cosmic variance of the kurtosis field. 



4 Discussion 

In the present paper we showed how to implement the multiple impulse approximation for 
perturbations on a photon beam (stemming from the effect of the string network) in the 
actual construction of higher order correlations for the CMB anisotropics. We then focused 
on the four-point function and on the excess kurtosis parameter, finding for the latter a value 
K. ~ 10- 2 . 

We also calculated explicitly the rms excess kurtosis ICcv predicted to exist even for a 
Gaussian underlying field and showed its dependence on the primordial spectral index of 
density fluctuations. This constitutes the main source of theoretical uncertainty at COBE 
scales. In fact, the cosmic string signature that might have been observable is actually 
blurred in the cosmic variance mist reigning at very large scales. 

Nevertheless, there is still a chance of getting a string-characteristic angular dependence 
from the study of the mean four-point correlation function by exploiting the particular 
geometries deriving from it; namely, its collapsed cases (where some of the five independent 
angles are taken to be zero) or from particular choices for these angles (as in the case of 
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Figure 1: Excess Kurtosis parameter of a Gaussian temperature fluctuation field as function 
of the spectral index. 

taking all angles equal). 

A preliminary analysis (Gangui & Perivolaropoulos 1994) making use of just one non- 
vanishing angle in a collapsed configuration as the one mentioned above shows a potentially 
interesting effect that could eventually increase notably the small non-Gaussian signal, and 
suggests that this is indeed a subject worth of further investigation. Some of these alterna- 
tives are presently under study and we expect to report progress on this subject in a future 
publication. 
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